Classifying Novel Phases of Spinor Atoms 
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We consider many-body states of bosonic spinor atoms which, at the mean-field level, can be 
characterized by a single-particle wave function for the BEC and Mott insulating states. We describe 
and apply a classification scheme that makes explicit spin symmetries of such states and enables 
one to naturally analyze their collective modes and topological excitations. Quite generally, the 
method allows classification of a spin F system as a polyhedron with 2F vertices. We apply the 
method to the many-body states of bosons with spins two and three. For spin-two atoms we find 
the ferromagnetic state, a continuum of nematic states, and a state having the symmetry of the 
point group of the regular tetrahedron. For spin- three atoms we obtain similar ferromagnetic and 
nematic phases as well as states having symmetries of various types of polyhedra with six vertices: 
the hexagon, the pyramid with pentagonal base, the prism, and the octahedron. 
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Ultracold atoms in either a single optical trap or in 
an optical lattice provide clean realizations of unique 
systems of spins which were previously studied only as 
toy mathematical problems (for a recent review, see 
(H). Depending on which hyperfine state is popu- 
lated, alkali atoms can have spin one or two. Differ- 
ent phases of spin-one alkali bosons have been experi- 
mentally realized 0, 0, 0, 0] and considered theoretically 
S S B uL OH El El • Spin- two bosons have also been ex- 
perimentally probed [1, H, [H, [3] and theoretically stud- 
ied for the case of a single optical tra p ITR [l6[ as well 
as, very recently, an optical lattice [13, Il8|. Finally, the 
Stuttgart group has recently succeeded in obtaining a 
Bose-Einstein condensation of 52 Cr atoms [HI which are 
spin-three bosons. This was followed by theoretical work 
13, Hi] showing the possible types of phases that can be 
realized for such a spin-three system. 

Classification schemes of single-particle states with 
non-zero spins are needed to describe both superfluid 
condensates and Mott insulating states of spinor bosonic 
atoms. However, such classification becomes increasingly 
more difficult for larger spins. For instance, classifying 
the state of a spin-half particle is straightforward; only 
knowledge of the expectation value of the spin operator 
(F) is needed. On the other hand, for a spin-one particle, 
knowledge of the expectation value of the nematic tensor 
familiar from the classical theory of liquid crystals [22| 
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Qab = -^{F a F h + F b F a ) — -5 a bF 
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in addition to (F) is required. Proceeding along these 
lines, one finds that for larger spin such a classification 
scheme becomes quite cumbersome since one needs to 
consider order parameters that involve higher-order prod- 
ucts of spin operators, and a physical interpretation is not 
immediate. 

In this Letter, we present an alternative classification 
scheme which will work well for large spin. This scheme 
allows the symmetries of a general spin F particle to be 



represented by a polyhedron with 2F vertices. To illus- 
trate the method we use it to discuss spin-two bosons 
which are naturally realized as a hyperfine state of alkali 
atoms. The previously discussed superfluid ferromag- 
netic, polar, and cyclic condensates [15|, 16] are shown 
to have the transformation properties of a ferromagnet, 
a nematic (either uniaxial or biaxial), and a tetrahedron. 
In addition, we show that the nematic phase has an addi- 
tional degeneracy at the mean-field level. We then show 
that states having similar transformation properties will 
be obtained in the Mott insulating state. Next, we show 
that the classification scheme immediately gives the num- 
ber of Goldstone modes and allows one to classify topo- 
logical excitations. In fact, in the biaxial nematic phase 
and the tetrahedratic phase, the topological excitations 
combine according to nonabelian groups [23]. Finally, 
due to recent experimental interest, we will also discuss 
classification of spin-three bosons. 

We will now outline our general classification scheme 
for a single spin. The key ingredient of our classification 
scheme is identifying states of spin F particles with 2F 
points on the unit sphere j24[. Consider a particle of 
spin F in the state given by = J2 a= _ F A a \a) where 
F z \a) = a\a) and A a are a set of normalized complex 
coefficients. To gain a physical understanding of this 
state, the idea is to find the set of "maximally polar- 
ized" states which are orthogonal to The maximally 
polarized state \() pointing in the n = (0,<f) direction 
is determined by the equation F • h\Q = F\(). A conve- 
nient (not normalized) representation of this state (which 
is related to the Schwinger boson representation) is 
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a=0 



where ( = e^tan(#/2) is the stereographic mapping of 
the unit sphere to the complex plane. We define the 



characteristic polynomial for in £ to be 



MO 



(3) 



The values of the 2F complex roots of /(C), {0} = {C : 
/0(C) = 0}> (which are in one-to-one correspondence with 
a set of points on the unit sphere {(6^,<^)}) determine 
the coefficients A a and therefore up to normalization 
and an overall phase factor. Most importantly, the sym- 
metries of correspond to the operations under which 
the set of points on the unit sphere {(^,^)} are invari- 
ant. A similar method has been developed in the 19th 
century mathematics community to solve quintic polyno- 
mials in terms of rotation of regular icosahedra [25[ . 

We illustrate our general approach by considering a 
specific problem of spin- two bosons, first concentrating 
on the BEC state in a single optical trap and then ex- 
tending the discussion to the insulating Mott state in an 
optical lattice. Spin two bosons interact with the contact 
potential 

y int ( Xl - x 2 ) = <S(xi - x 2 )(^ ^o + 92V2 + £4^4) (4) 

where Vf projects into the state with total spin F and 
gp = l-Kaplvn where ap is the scattering length cor- 
responding to spin F. Note the requirement that the 
wave functions be symmetric under interchange of parti- 
cles prohibits odd spin projection operators. In the su- 
perfluid phase when all atoms occupy the same orbital 
state, potential (|4j) leads to the interaction Hamiltonian 
of 



state, poterj 
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Hint = ^U n(n 



1) + \uiVq + ^U 2 (F 2 



6n) (5) 



where Uo = y(4<22 + 3*24), U± = y(7ao — 10a2 + 3*24), and 
U2 = ^(04 — 02) (7 is a positive constant). In Eq. ([5]) 
n = o) a a a (here and after sums over repeated indices 
are implied) where a^ a creates a spin-two boson in an 
eigenstate of F z having eigenvalue a. Finally, F 2 = F • F 
where F = a^Tapap and T a p are the spin- two matrices. 

The mean-field description for a BEC of TV bosons is 
|^Sf) = ^=y (^ofla) |0), where the normalized com- 
plex coefficients A a are the variational parameters. Since 
the kinetic energy of the system will not depend on 
the spin degrees of freedom, it is sufficient to minimize 
(V ; SF|7^int|V ; SF) over the variational parameters. To lead- 
ing order in TV, we find with this wave function 



(H int )/N 2 = — U 1 \2A 2 A_ 2 -2A 1 A_ 1 + A A \ 2 
+ lu + lu 2 (A* a T a pA ) 2 . 



(6) 



Carrying out the minimization with the constraint of nor- 
malization, A^A a = 1, we find the phase diagram given 
in Fig. [U where the boundaries are given by 02 — cla = 



Ferromagnetic 
(1,0,0,0,0) 




FIG. 1: Phase diagram for spin- two bosons in a single optical 
trap. The axes are in units of as- For the Mott insulating 
state with one bosons per site, the phase diagram will be the 
same with horizontal and vertical axes corresponding to eo — £4 
and 62 — £4 respectively (see Eq. 112)) . 



±^(ao — cla) and 02 =0,4. The phases correspond to the 
coefficients (A_2, ^--1, A)? A.2) given by 

sin(?7) 



T: (1,0,0,0,0) 
'sin(ry) 




^_ ,0,cos(t7),0, ^_ 



(7) 
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which are uniquely determined up to SO (3) spin rota- 
tions by the minimization. Phase M has the additional 
degeneracy parameter n which does not correspond to 
a symmetry of the hamiltonian. One should note, how- 
ever, that this degeneracy only arises at the mean-field 
level and will be removed when we include corrections to 
the mean-field wave function j26[. On the other hand, in 
real experiments, the quadratic Zeeman effect provides 
another way of removing this degeneracy which we will 
discuss. Using the scattering lengths taken from Ref. 
[lH we also indicate where 23 Na and 87 Rb will lie in the 
phase diagram. The same phase boundaries were found 



by Ciobanu, Yip, and Ho 



15[ | for the superfluid phase. 



However, we emphasize that there is a continuous degen- 
eracy in the phase J\f at the mean-field level. 

Symmetries of the phases given by the wave functions 
in Eq. Q is not immediate upon first glance. We will 
now illustrate the utility of our method by classifying the 
various phases. For spin two, the maximally polarized 
state as in Eq. ([2]) is 

|C) = C 4 |"2) + 2C 3 |-1) + V6C 2 |0) + 2C|1) + \2). (8) 

Phase T is clearly the ferromagnetic state. However, it 
will prove instructive to proceed with our systematic clas- 
sification scheme. The characteristic polynomial Eq. ([3]) 
for this case is simply f^ T (£) = (" 4 which has the four- fold 
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degenerate root of ( = which corresponds to the north 
pole of the unit sphere (0, (j)) = (0,0). The only symme- 
tries this state possesses are therefore the U(l) group of 
continuous rotations about the z-axis. 

We now move on to phase J\f. Before applying our clas- 
sification method, it is first useful to consider the order 
parameters constructed in the conventional way. First, 
we note that (F) = 0. Next the eigenvalues of the ne- 
matic order parameter Eq. (pQ) are 



-2cos(2t?) 
Eig(Q a6 ) = ( cos(277) + >/3sin(277) 
cos(2t7) — v // 3sin(2^) 



(9) 



Now, applying our method the characteristic polynomial 
is 

UAO = \j\ sm(r?)C 4 + VScosfaK 2 + \J\ sinfa). (10) 

The roots of this equation all lie at the vertices of a 
rectangle. For < 77 < 7r/3, tt/3 < 77 < 27r/3, and 
27r/3 < 77 < 7r the rectangle will lie in the zy, xy, and 
xz planes respectively. We find that 77 = nix) 3 (defined 
ATI) corresponds to two sets of double roots at opposite 
poles. This state will have the transformation properties 
of a uniaxial nematic. That is, there is the symmetry of 
continuous rotation about the nematic axis. There are 
also the additional symmetries of discrete rotations by tt 
about any axis in the plane perpendicular to the nematic 
axis. On the other hand, when 77 = (n + l/2)7r/3 (defined 
Af2) the roots lie at the vertices of a square which has no 
continuous symmetries (contrary to what one might be 
led to believe from Q a b)- The symmetries here are rota- 
tions by 7r/2 about the nematic axis and rotations by tt 
about the other two principal axes. For all other values of 
77 the nematic will be biaxial, which we refer to as phase 
Af3. This case has the fewest symmetry operations, the 
rotations by tt about the three principal axes. 

A route to remove the mean-field degeneracy associ- 
ated with 77 is to apply an external magnetic field. Since 
for typical experiments the time it takes a spin state to 
relax is longer than the trap lifetime Q , it is appropriate 
to think of (F z ) as a conserved quantity, and the linear 
Zeeman term is unimportant. Thus, the most impor- 
tant effect due to an external magnetic field is due to 
the quadratic Zeeman term (proopertional to F%) which 
can either have a positive or negative coefficient. In the 
(most common) event that the coefficient is negative, the 
phase Af2 which transforms as the square will be favored. 
On the other hand, if the coefficient is positive, the Afl 
phase will be favored. 

We now consider the remaining region in the phase di- 
agram, phase T, which, as we will see, is tetrahedrat- 
ically ordered. For this case, Q a b gives no useful in- 
formation since (Q a b) = as well as (F) = 0. Thus, 
the utility of our classification method will be demon- 
strated most clearly here. The characteristic polynomial 



for this situation is fip T (() = \J \C A + 2y §C The roots 
of this equation correspond to the points {(^,^)} = 
{(O,O),(0,tt/3),(0, tt),(<9, 5tt/3)} where = 2tan- 1 (v / 2). 
These points lie at the vertices of a regular tetrahedron 
on the unit sphere. Thus \i/jt) nas the full symmetry of 
the point group of the tetrahedron. 

We will now consider what happens when the spin- 
two bosons are in an optical lattice in a Mott Insu- 
lating state (one per site). Such a system will be de- 
scribed by the hamiltonian TL = 7Ykin + T^Hub where 
Hkin = -J^2(ij)( a la a j^ + h.c.) describes hopping be- 
tween adjacent sites and Hiiub gives the on-site interac- 
tion 



^Hub,i = -U ni(ni-1)- 



lu.V^MF?-^). (11) 



Note the similarity between the on-site interaction and 
Eq. ([5]). When the sites are completely decoupled, the 
ground state energy will be independent of the spin state 
of any given site leading to a macroscopic degeneracy 
which will be removed by hopping. Treating the hopping 
perturbatively, we arrive at the effective hamiltonian 



H eff = ( e oPo(ij) + e 2 V 2 (ij) + e A V A (ij)) 



(12) 



C/0+E/1-6C/2 ' 62 " u -3U 2 ' 
This effective hamiltonian is valid in 



which acts in the subspace of singly occupied sites. Here, 
Vf{^3) projects neighboring sites at ij into a state with 
total spin total spin F and eo _ 

and e 4 = Uo +i U2 
the limit J <C Uq. We point out that 7i e fi can also be 
written in terms of powers up to four of spin operators 
on adjacent sites • Sj, but we find the representation 
in Eq. ([T2]) easier to work with. 

When all of the scattering lengths are equal, we will 
have for the on-site parameters U\ = U2 = 0, which will 
make eo = 62 = 64. This corresponds to the situation 
where the hamiltonian ([T2]) has a special SU(5) symme- 
try. For this situation, the state with all sites having the 
same spin wave function I^mi) = n^a|<^)i will be an 
exact eigenstate of H e ff, which is the ground state. There 
will still be a rather large degeneracy for this situation 
since the ground state energy will not depend on the co- 
efficients A a . On the other hand, when the scattering 
lengths begin to differ this degeneracy will be partially 
lifted. We will now determine the phase diagram using 
the variational wave function. By evaluating (^l^effl^) 
and minimizing over the set of five complex coefficients 
A a we find precisely the same coefficients as in flJJ) . More- 
over, the phase diagram for this case will be identical 
with that given in Fig. [T] with the the axes changed to 

At the level of virtual hopping Eq. ([T2]) in perturba- 
tion theory (second order in J/U), the states in phase M 
are degenerate with respect to the parameter 77. How- 
ever, this does not correspond to a true symmetry of the 



parent hamiltonian. One therefore expects that this de- 
generacy will be removed at the ring exchange level in 
perturbation theory (fourth order in J/U). On the other 
hand, quantum fluctuations at the virtual hopping level 
from 7Y e ff may favor particular phases which can be of 
the same order as ring exchange. The detailed treatment 
of this problem is beyond the scope of the present work 
and will be considered elsewhere [261 ]. 

Knowing the symmetry properties of the wave func- 
tion leads to a natural classification of the collective ex- 
citations of the system. Breaking continuous symmetries 
leads to Goldstone bosons; the number of such modes 
corresponds to the number of broken symmetry genera- 
tors. Our hamiltonians will always have an SO (3) sym- 
metry which corresponds to global spin rotations (at the 
point where all the scattering lengths are the same there 
will be a special SU(5) symmetry). As can be seen in 
the pictorial representation, the uniaxial nematic phase 
Ml has one remaining continuous symmetry and thus 
break two of the generators of SO (3). This state will 
therefore have two Goldstone bosons. The biaxial ne- 
matic phases M2 and A/*3, and tetrahedratic phase T 
all have no remaining continuous symmetries, so these 
will have three Goldstone modes. The topological excita- 
tions out of the various phases can also be classified using 
the homotopy theory. The fundamental group describes 
how different defects combine. For instance, it is known 
that the fundamental group of the biaxial nematic J\f3 is 
the (nonabelian) quaternion group 0. The phase M2 
(transforming as the square) has the larger 16-element 
dicyclic group J2?| for its fundamental group, giving it 
more types of topological excitations. Finally, the fun- 
damental group of the tetrahedratic state is the binary 
tetrahedral group with 24 elements [27[. The number of 
topological excitations (4, 6, and 6 respectively) is the 
number of nontrivial conjugacy classes in each of these 
groups. 

We now will briefly comment on the types of states 
one can expect for the spin-three problem motivated by 
recent interest [l9|, [lo|, A similar variational wave 

function approach for either the superfluid or Mott insu- 
lating states as done above can be used for the spin-three 
case. The characteristic polynomial /^(C) will have six 
roots (instead of four as in the spin- two case). Possible 
phases that can be realized for such a system are shown 
in Fig. [2j The boundaries between these phases can be 
found in H3, EH; the purpose here is to show that our 



classification scheme can be applied naturally to the spin- 
three system. 

In conclusion, we have presented a classification 
scheme for cold spinor bosons. We considered spin- two 
bosons as an illustrative example, which can be realized 
as a hyperfine state of an alkali atom. In addition to 
the ferromagnetic and nematic states, we showed that a 
tetrahedratic ordering can occur for this system. Finally, 
we discussed an extension to spin-three bosons. 
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FIG. 2: Some possible phases that can be realized for a 
spin-three system either in the superfluid or Mott insulat- 
ing states. Shown are the coefficients of the wave functions 
A a and the shapes representing the symmetries of the wave 
functions. The phases transform as the following polyhedra 
(a) the hexagon, (b) the pyramid with pentagonal base, (c) 
the prism, and (d) the octahedron. For phase (c) we have 
the condition tan 2 (77) < 10. Additional phases similar to the 
spin two case (not shown) are the ferromagnetic and nematic 
states. 
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